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Abstract. The interpretation of some electrovacuum spacetimes in terms of 
counterrotating perfect fluid discs is presented. The interpretation is mades by means 
of an "inverse problem" approach used to obtain disc sources of known static solutions 
of the Einstein-Maxwell equations. In order to do such interpretation, a detailed 
study is presented of the counterrotating model (CRM) for generic electrovacuum 
static axially symmetric relativistic thin discs with nonzero radial pressure. Four 
^ • simple families of models of counterrotating charged discs based on Chazy-Curzon-type, 

bX), Zipoy-Voorhees-type, Bonnor-Sackfield-type, and charged and magnetized Darmois 

electrovacuum metrics are considered where we obtain some discs with a CRM well 
behaved. 
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1. Introduction 

X 

■ Several methods are known to exactly solve the Einstein and Einstein-Maxwell field 
equations, or to generate new exact solutions from simple known solutions |1] . However, 
the above mentioned methods in general lead to solutions without a clear physical 
interpretation or to solutions that depend of many parameters without a clear physical 
meaning. Accordingly, it is of importance to have some appropriate procedures to obtain 
physical interpretations of these exact solutions. So, in the past years such procedures 
have been developed for static and stationary axially symmetric solutions in terms of 
thin and, more recently, thick disc models. 

Stationary or static axially symmetric exact solutions of Einstein equations 
describing relativistic thin discs are of great astrophysical importance since they can 
be used as models of certain stars, galaxies and accretion discs. These were first studied 
by Bonnor and Sackfield [2] , obtaining pressureless static discs, and then by Morgan and 
Morgan, obtaining static discs with and without radial pressure [31 H] . In connection 
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with gravitational collapse, discs were first studied by Chamorro, Gregory, and Stewart 
[5]. Discs with radial tension have been also studied [6]. 

Several classes of exact solutions of the Einstein field equations corresponding 
to static and stationary thin discs have been obtained by different authors, with or 
without radial pressure (see refs. to [IZj). Except for the pressureless disks all 
the other disks have as source matter with azimuthal pressure (tension) different from 
the radial pressure (tension). However, in some cases these disks can be interpreted 
as the superposition of two counterrotating perfect fluids. A detailed study of the 
counterrotating model for the case of static thin disks is presented in [I8|. Recently, 
more realistic models of thin disks and thin disks with halos made of perfect fluids were 
considered in [T9] . 

Disc sources for stationary axially symmetric spacetimes with magnetic fields are 
also of astrophysical importance mainly in the study of neutron stars, white dwarfs and 
galaxy formation. Although discs with electric fields do not have clear astrophysical 
importance, their study may be of interest in the context of exact solutions. Thin discs 
have been discussed as sources for Kerr-Newman fields pO], magnetostatic axisymmetric 
fields |2T], conformastationary metrics [22], and electrovacuum static counterrotating 
discs [23]. Following Ref. [20] the resulting discs can be interpreted either as rings with 
internal pressure and currents or as two counterrotating streams of freely moving charged 
particles, i.e., which move along electrogeodesics (solution to the geodesic equation in 
the presence of a Lorentz force). 

In all the above cases, the discs are obtained by an "inverse problem" approach, 
called by Synge the '^g-method^ [21]. The method works as follows: a solution of the 
vacuum Einstein equations is taken, such that there is a discontinuity in the derivatives of 
the metric tensor on the plane of the disc, and the energy-momentum tensor is obtained 
from the Einstein equations. The physical properties of the matter distribution are 
then studied by an analysis of the surface energy-momentum tensor so obtained. Very 
recently the above procedure has been generalized in order to obtain thick disc models 
[25] and also non-axisymmetric planar distributions of charged dust [26]. As we can see, 
by this "inverse problem" approach we can do physical interpretations of many known 
solutions of the Einstein and Einstein-Maxwell vacuum equations, in the sense that the 
thin disks can act as exact sources for the space-time metrics given by the vacuum 
solutions. 

The purpose of this paper is the interpretation of some static electrovacuum 
spacetimes in terms of counterrotating perfect fluid discs. In order to do such 
interpretation, a detailed study is presented of the counterrotating model (CRM) for 
generic electrovacuum static axially symmetric relativistic thin discs with nonzero radial 
pressure. The counterrotating model for the case of static thin discs with radial pressure 
was recently studied in [18j and for static discs without radial pressure but without 
electromagnetic fleld in [23], so the material presented here is a continuation of the 
mentioned works. The paper is structured as follows. In Sec. |2], we present a summary 
of the procedure to construct models of discs with nonzero radial pressure. We obtain 
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expressions for the surface energy-momentum tensor and the surface current density of 
the disc. 

In Sec. El the discs are interpreted in terms of the CRM. We find a general 
constraint over the counterrotating tangential velocities needed to cast the surface 
energy-momentum tensor of the disc as the superposition of two counterrotating charged 
perfect fluids. We also find expressions for the energy density, pressure, current density, 
and tangential velocity of the counterrotating fluids. We show that, in general, it is not 
possible to take the two counterrotating fluids as circulating along electrogeodesics nor 
take the two counterrotating tangential velocities as equal and opposite. 

In Sec. m four simple families of models of counterrotating charged discs based 
on Chazy-Curzon-type, Zipoy-Voorhees-type, Bonnor-Sackfield-type, and charged and 
magnetized Darmois metrics are presented where we obtain some discs with a well 
behaved CRM. In particular, we study the tangential velocities, energy and electric 
charge densities in the particular case when both streams have equal pressure or equal 
energy density. Also the stability against radial perturbation is analyzed. Finally, in 
Sec. [5|, we summarize our main results. 



2. Electrovacuum static relativistic discs 

2.1. Models of discs with nonzero radial pressure 

In this section we present a summary of the procedure to obtain electrovacuum static 
axially symmetric relativistic thin discs with nonzero radial pressure. The line element 
of a static axisymmetric field can be written in quasi-cylindrical coordinates {t,ip, r, z) 
in the form 

where A and R are functions of r and z only. The vacuum Einstein-Maxwell 
equations, in geometrized units in which SvrG = c = /iq = = 1; are given by 

Rah = Tah, (2) 

i^"';fe = 0, (3) 

with 

Tab = FacF^^ — -QabFcdF^'^, (4) 

Fab = daA, - d,Aa. (5) 

In order to obtain a solution of ([2]) - ([3]) representing a thin disc at 2 = 0, we 
assume that the components of the metric tensor are continuous across the disc, but 
with first derivatives discontinuous on the plane z = 0, with discontinuity functions 



Jab 



9ab,z\ r,+ 9ab,z\ f.- 2 Qab,! 



Thus, the Einstein-Maxwell equations yield an energy-momentum tensor Tab = ^at™ + 
TX\ where = Qab S{z), and a current density Ja = Ja5{z) = -2e2(*-^)Aa,^5(z), 
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where 6{z) is the usual Dirac function with support on the disk. T^^b" 
electromagnetic tensor defined in Eq. (jl]), ja is the current density on the plane z = 0, 
and 

Qt = i^ih^'K - h^'^^t + rhl - g'^K + Kia^'^t - r^l)} 

is the distributional energy-momentum tensor. 

The "true" surface energy-momentum tensor (SEMT) of the disc, Sab-, and the 
"true" surface current density, j^, can be obtained through the relations 



Sab 
Ja 



rpms.t J 
-I- ab 

J n dSri^ 



Q 



ab 



e^-" Ja 



(6) 
(7) 



where dsn = \Jgzz dz is the "physical measure" of length in the direction normal to the 
disc. For the metric ([1]), the nonzero components of S*^ are 



cO 

Dq — 


2e* 


-A 


sl = 


2e* 


-A 


sl = 


2e* 


-A 



A - 2^ + 

*-)2 ^ ^ 52 



A 



1Z 1 

R,z 



(8) 
(9) 
(10) 



and the components of the surface current density are equal to 

j, = -26-^-^,,, (11) 

-2e'^-M„ (12) 

where ip and A are the electrostatic and magnetostatic potentials, respectively, which 
are also functions of r and z only. All the quantities are evaluated at z = O"*". 

In order to give physical significance to the components of the energy-momentum 
tensor we project it onto the orthonormal tetrad ea'' = {V^^, W^, X^, Y^}, where 



= e-*(l, 0, 0, 0), X" = e^-^{0, 0, 1, 0), 



(13) 



1^^ = — (0,1,0,0), 



ya ^ e*-^(0, 0,0,1). 



In terms of this tetrad (or observer ) the metric, the SEMT and the current density may 
be decomposed as 



dab 
Sab 
Ja 



- VaVh + WaWh + XaX, + r.n, 

eVaVb + P^WaWb + PrXaXb, 
aVa+jWa, 



where 



Sq, Pip 



Si 



are, respectively, the energy density, the azimuthal pressure, the radial pressure. 



(14) 
(15) 
(16) 

(17) 

(18) 
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are the electric charge density and the azimutal current density of the disc measured by 
this observer with four-velocity V"', and 

g = e + + Pr (19) 

is the "effective Newtonian density" . 

2.2. Solution of the field equations 

For the metric ([I]), the Einstein- Maxwell equations in vacuum imply that R satisfies the 
Laplace's equation 

R ~\~ R zz ~ 0, (20) 

so that the function R can be considered as the real part of an analytical function 
F[iy) = R{r, z) +iZ{r, z), where u = r + iz. Thus the function -F(z/) defines a conformal 
transformation 

r^R{r,z), z^Z{r,z), (21) 

in such a way that the metric ([1]) takes the Weyl form 

ds^ = - e^^dt" + e-'^^fi^'V + e'\rfi?2 + dZ% (22) 

where R and Z are the Weyl's canonical coordinates. 

In these coordinates, the field equations ([2]) - (j3]) are equivalent to the usual complex 
Ernst equations [27J 

fA£ = {VS + 2<I>*V$) ■ VS, (23) 
/A$ = (V^ + 2<I>*V$) ■ V$, (24) 

where A and V are the standard differential operators in cylindrical coordinates, / = e^'^ 
and S = S* (static spacetime). The metric functions are obtained via 

/ + (25) 
" 1 
.4/' 

where V2( = R + iZ, so that -\/29,j = d^R — id^z- The electromagnetic potentials are 
related to $ via 

\X = v^^|(Im$),c, (27) 
and At = V2Re^. 

It is known that in static electromagnetic fields there is a lineal relation between 
the electrostatic potential and its magnetic counterpart [28], |29] . Some authors therefore 
study only electric fields. In this paper we consider the field equations when both electric 
and magnetic fields are presents, since the true magnetic potential will be different, and 
therefore the physical properties will differ too. 



:(£,f + 2<l>*$,(^)(^,^ + 2$$*^) - 



(26) 
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Once a solution of the system of equations (1251) - is known, we can obtain a 
solution of the field equations ([2]) - ([3]) in the original coordinates by setting 

(28) 
(29) 
(30) 
(31) 
(32) 

where F' = dF/ dv. With this transformation can be written as 

Vr = I =-FT- (33) 

IRI + R\ R 



R{r, 


z) 


= ReF{u), 


\t'(r, 


^) 


= f^{R,z), 


A(r, 




= X{R,Z) + \n\F'{iy 


ilj{r, 




= A{R,Z), 


A{r, 


z) 


= A^iR,Z), 



Thus, the construction of thin discs models with nonzero radial pressure depends 
strongly of the analytical function Fiy), due to the fact that this function gives the 
behavior of the radial pressure and also defines the conformal transformation that leads 
the field equations to the form given for Ernst (123|) - ( 12^ . A common choice for F{v) 
is 

F{u) = u + izo , (34) 

where zq is a constant. In this case we have infinite discs with zero radial pressure. This 
is the well known "displace, cut and reflect" method used in almost all of the works 
about discs models to generate thin discs from solutions of the Einstein and Einstein- 
Maxwell equations. Another possible choice for -F(z^) was presented in reference [6] and 
is given by 

F{u) = V + a - 1 , (35) 

where a > 0. With this choice for Fiv^ the general expresion for the radial pressure 
leads to 

_ 2ae^-^ 

[l + («2_ 1)^2] 1/2' ^-^^^ 

so that we obtain thin discs with nonzero radial pressure and of finite radius, located at 
2r = 0, 0<r<l. In order to obtain discs with non-unit radius, we only need to make 
the transformation r — > ar, where a is the radius of the disc. We shall also take F{y) 
as given by ( l35l) . 

With this choice for F(v) the image of the disc by the conformal mapping (12 ip is 
the surface 

c?R^ ^Z'' = c? (37) 

so that the discs are mapped into spheroidal thin shells of matter and its exterior is 
mapped into the exterior of the shells. We have three possible values for a: a = 1, a 
spherical shell, a > 1, a prolate spheroidal shell, and < a < 1, an oblate spheroidal 
shell. By considering the above shells as sources, we then seek for exterior solutions of 
the field equations (l23l) - (!24l) adapted to the symmetry of the shells. Then, using (!28ll 
- (l32l) we obtain the corresponding disc solutions in the original coordinates. 
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3. The counterrotating model 

3.1. Counterrotating charged perfect fluid discs 

We now consider, based on Refs. [30] and [31], ttie possibility that the SEMT 5*"^ and 
the current density j"" can be written as the superposition of two counterrotating charged 
perfect fluids that circulate in opposite directions; that is, we assume 

gab ^ gab ^ gab ^ (gg^, 
/ =j++j-, (39) 

where the quantities on the right-hand side are, respectively, the SEMT and the current 
density of the prograde and retrograde counterrotating fluids. 

Let = (f/° , f/j_, 0, 0) be the velocity vectors of the two counterrotating fluids. 
In order to do the decomposition (|38l) and ( |39l) we project the velocity vectors onto the 
tetrad ea^, using the relations [32] 

f/| = e\Ui, Ul = Bf^^ui. (40) 

In terms of the tetrad (fT3l) we can write 



Ul = I , . (41) 
'1 — vi 



where v± = U^/U^ are the tangential velocities of the streams with respect to the 
tetrad. 

Another quantity related with the counterrotating motion is the specific angular 
momentum of a particle rotating at a radius r, defined as h± = g^^U±. We can write 

= (42) 



This quantity can be used to analyze the stability of the discs against radial 
perturbations. The condition of stability, 

^ > 0. (43) 

dr 

is an extension of Rayleigh criteria of stability of a fluid in rest in a gravitational field 



In terms of the metric hab = Qab — YaYb of the hypersurface 2; = 0, we can write the 
SEMT as 

5"^' = {e + pr)V''V'' + (p^-p.)iyW' + (44) 
and using (14T!) we obtain 

^,,^f_(^.,v.){l-vl)UlUX 



+ 



(f+ — f_)^ 

f{v+,v+){l-vl) U^Ut 
(f+ — f_)2 
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+ Prh 



ab 



where 



f{Vi,V2) = {e + Pr)ViV2 + P^ - Pr- (45) 

Clearly, in order to cast the SEMT in the form fl38|) . the mixed term must be absent 
and therefore the counterrotating tangential velocities must be related by 

= 0, (46) 

where we assume that \v±\ 7^ 1. 

Then, assuming a given choice for the tangential velocities in agreement with the 
above relation, we can write the SEMT as fl38|) with 

St = {e±+p±)UiUi + p±h-\ (47) 

so that we have two counterrotating perfect fluids with energy densities and pressures, 
measured by an observer comoving with the streams, given by 



e± + P± 
P++ p_ 



{e + p,)v^, (48) 

Pr. (49) 

From these equations also follows that 

e+ + e_ = e + Pr - p^- (50) 

Thus the SEMT 5*"* can be written as the superposition of two counterrotating perfect 
fluids if, and only if, the constraint ( 146|) admits a solution such that f+ 7^ v^. This 
result is completely equivalent to the necessary and sufficient condition obtained in Ref. 

Similarly, substituting fj4T|) in f|T6|) we can write the current density as fl39l) with 

j± = <y±ui, (51) 

where a± are the counterrotating electric charge densities measured by an observer 
comoving with the streams, which are given by 



1 



4 



(J 



(52) 



Thus, we have a disc makes of two counterrotating charged perfect fluids with energy 
densities and pressures given by and (j^ . and electric charge densities given by 
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3.2. Counterrotating tangential velocities 

All the main quantities associated with the CRM depend on the counterrotating 
tangential velocities v±. However, even for known values of v±, the system of equations 
f HHj) - (H9l) is underdetermined so that we can have information about the physical 
quantities corresponding to each stream only in some particular situations. We will 
examine two simple cases in this paper, when the two streams have equal pressure and 
when they have equal energy density. On the other hand, the constraint ( l46l) does not 
determine v± uniquely so that we need to imposse some additional requeriment in order 
to obtain a complete determination of the tangential velocities leading to a well defined 
CRM. 

A possibility, commonly assumed, is to take the two counterrotating streams as 
circulating along electrogeodesics 
1 



-e±gab,rUiU^ = -a^FraU^. 



(53) 



Let uj± = U^/Uj_ be the angular velocities of the particles. In terms of u^. fl53l) takes 
the form 



1 / 2 

-e±{gu^r^^± + goo,, 



and v± can be written as 



1 



UJ±. 



(54) 



(55) 



As an example, let us consider the electrogeodesic motion of the particles when 
we have equal energy density. By simplicity, we will analyze the electrostatic and 
magnetostatic cases separately. The angular velocity uj± can be obtained from the 
electrogeodesic equation (IMI) and, in the electrostatic case, we obtain u± = ±u with 



,4* 



R 



-2^,,(^,, -A,^)+V',rV',.e 



-2'I' 



-2(^,^ - A^^){R^r - R'^^r) + Rlp,r'ip,ze~ 



2* 



and therefore 



-R 



2(\1' 2 — A^z)dpr/dr + Prljj^r^^,^ 



-21' 



(56) 



(57) 



-2(^,^ - A,,){R,r - R^,r) + Rlp,r'ip,ze-'^^ ^ 

where we have used the Einstein-Maxwell equations ([2]) and ([3]) for A and the expressions 
dH]) - ffT^ for the SEMT and the current density. It follows immediately that /(w+,f_) 
vanishes if 



-2il' 



1 dpr 

Jr'^^ 2(^,,-A,,)' 



and therefore 



Ji{r) 



Pr = Cie" 
where Ci is a constant and 

fi{r) = - 



-2* 



2(^,.-A, 



-dr. 



(58) 



(59) 



(60) 



Counterrotating perfect fluid discs as sources of electrovacuum static spacetimes 10 



Similarly, for magnetostatic fields we obtain again uj± = ±00 with 



A.,) + A,rA„ 



,2* 



2i?(^,, - A,,)iRr - R'^,r) + A,rA,,e- 



21' 



so that 



R 



2* 



2(^,, - A^,)dpr/dr -prA,rA^,e 
-2(^,, - A,,){R,r - R^,r) + A^rA,,e^^/R 



That is, vanishes if 

1 dpr 



A rA 



2* 



Pr dr 



2^.- A 



and therefore 



J2{r) 



Pr = C2& 

where again C2 is a constant and 



/2(r) 



2* 



-(ir. 



(61) 



(62) 



(63) 



(64) 



(65) 



2(*.-A,., 

That is, for electrostatic or magnetostatic fields and equal energy densities the streams 
circulate along electrogeodesics only if the radial pressure satisfies f l59|) or 
respectively. Note that the two electrogeodesic fluids are circulating with equal and 
opposite tangential velocities. In more general situations, for example with simultaneous 
presence of electric and magnetic fields, Pr is a complicated function of r that can be 
obtained by a similar procedure and the streams would move with different velocities. 

As we can see, in order to have the counterrotating fluids circulating along 
electrogeodesics, the radial pressure must satisfy a very strong condition that, in general, 
is not satisfied by the thin disc solutions. Accordingly, we need to consider an additional 
requeriment different of the electrogeodesic motion in order to exactly determine the 
counterrotating velocities. Another possibility is to take the particles of the streams 
not circulating along electrogeodesics but with equal and opposite tangential velocities, 
that is 



so that fj46|) is equivalent to 

-Pr 



e + pr 



(66) 



(67) 



This choice, commonly considered, leads to a complete determination of the velocity 
vectors. However, this can be made only when the above expression is positive definite. 
In particular, when the two streams have equal pressure or equal energy density we have 
p± = \pr and e± = |(e + pr - p^). 

In the general case, when (1671) is not positive definite, the two fluids circulate with 
different velocities and we can write fH6l) as 



Pr -Py 

e + pr 



(68) 
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and so we can obtain a CRM only if 

\Pr-Pv\ < \ e + pr \. (69) 

However, this relation does not determine completely the tangential velocities, and 
therefore the CRM is undetermined. 



4. Some simple examples of counterrotating charged discs 
4.I. Discs from a Chazy-Curzon-type solution 

The first family of solutions to consider is a Chazy-Curzon-type solution which is given 
by 

= 7 , , ^ r, (70) 



A^-^. (n, 



* {a + l)e"'/P-{a-l)e-"'/p' ^ ' 

= 72752 cos ^, (73) 

where = 1 -|- 6^, with 6^ = q\ + q^, and 7 is a real constant. Here qi and q2 are the 
electric and magnetic parameters, respectively, p and 9 are quasi-spherical coordinates, 
related to the Weyl coordinates by 

R = psm6, Z = pcosO, (74) 

where < p < 00 and — f < ^ < f . In this coordinates we must choose a = 1 
and the shell would be located at p = 1. This solution can be generated, in these 
coordinates, using the well-known complex potential formalism proposed by Ernst [27] 
from the Chazy-Curzon vacuum solution [311 ES], by choosing the parameter q of Ref. 
f27\ as complex. For ^2 = we have an electrostatic solution and for gi = one obtains 
its magnetostatic analogue [36] . 

From the above expressions we can compute the physical quantities associated with 
the discs. We obtain 

(75) 
(76) 
(77) 
(78) 



e 


= g-Pr-Pv 


Pv 


= Pr{l + 7^r^^ 


Pr 


= Poe^ ' , 




QiPo 


a 


V2 



where 



= (a + l)e.-(a-l)e-. ^''^ 



is a constant and 



6=^[ia + l)e^ + (a - l)e-^]pr. (80) 
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The strong energy condition requires that ^ > 0, so that we must choose 7 > 0. We can 
see that po is a positive constant for these values of 7. 

In order to study the behavior of these quantities we perform a graphical analysis 
of them for discs with 7 = 2 and different values of qi and q2, as functions of r. In Figs. 
[U^a) and [T](6) we show the energy density and the pressures pr and for qi = q2 = 0, 
0.5, 1.0, and 1.5. We see that the energy density initially is a positive quantity and 
then becomes negative in violation of the weak energy condition. Therefore only the 
central region of these discs has a physically reazonable behavior. We can observe that 
the pressures Pr and p^ are always positive quantities everywhere on discs. We also 
see that the presence of electromagnetic field diminishes the energy density near to the 
center of the discs and later increases it, and it diminishes the pressures everywhere on 
discs. Next, the charge density cr is represented in Fig. [l](c) for qi = q2 = 0, 0.2, 0.3, 
0.5, and 1.5. As we can see, it exhibits a similar behavior to the radial pressure. We 
also computed these functions for other values of the parameters (with 7 > 0) and, in 
all the cases, we found a similar behavior. 

We now consider the CRM. All the significant quantities can also be expressed in 
analytic form from the above expressions but the results are so cumbersome that better 
turns out only to analyze them graphically. Let us consider the case when the two fluids 
move with equal and opposite tangential velocities. In Fig. Wid) we plot the tangential 
velocity curves of the counterrotating streams f ^ for discs with 7 = 0.7 and qi = q2 = 2, 
3, 5, and 10. We see that the inclusion of electromagnetic field can make the velocities 
of the particles smaller than the light velocity. In Fig. [T](e) we have drawn the angular 
momentum h"^ for the same values of the parameters. In some cases we have strong 
changes in the slope at certain value of r, which means that there is a strong instability 
there. For other values of the parameters h'^ is an increasing monotonous function of r 
that corresponds to a stable CRM for the discs. Finally, in Fig. Wif) we have plotted 
the energy and electric charge densities e± and a± in the particular case when the two 
streams have equal pressure or equal energy density, for values of the parameters for 
which < 1. We see that a± behaves in the opposite way to a, whereas e± has a similar 
behavior to the energy density e. Therefore, for Chazy-Curzon-type fields we can build 
counterrotating thin disc sources only with a well-behaved central region. 

4-2. Discs from a Zipoy-Voorhees-type solution 

The second family of solutions considered is a Zipoy-Voorhees-type solution which can 
be written as 



e' 



2(x2 - 1)^/2 



(81) 



A 



(a + l)(x + l)T- (a - l)(x-l)T'' 

7^ 1 r — 1 

2 x'^ — y"^ 



(82) 



y2gi[(x + l)^-(a;-l)^] 



(83) 



(a + l){x + 1)T - (a - l)(x - 1)^' 
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= V2k-fq2y, (84) 

where = 1 + 6^, with = Q'l + and 7 is a real constant. Here qi and q2 are again 
the electric and magnetic parameters, respectively, x and y are the prolate spheroidal 
coordinates, related to the Weyl coordinates by 

= k'^{x^ -l){l-y^), Z = kxy, (85) 

where 1 < a; < 00 and < y < 1. In this case we must choose a > 1 and the shell 
would be located aX x = a/k > 1 and y = yl — with k = \/a^ — \. 

This solution can also be generated, in these coordinates, using the well-known 
complex potencial formalism proposed by Ernst [27] from the Zipoy-Voorhees vacuum 
solution [371 [38], also known as the Weyl 7-solution [391 HQ], by choosing the parameter 
q of Ref. [27] as complex. For g2 = we also have an electrostatic solution and for 
gi = one obtains its magnetostatic equivalent [36]. The case g2 = and 7 = 1 
corresponds to the Reissner- Nordstrom solution [T], in which case a = m/k, qi = e/k, 
with k"^ = m? — e^, so that = 1 + g^, m and e being the mass and charge parameters, 
respectively. When gi = and 7 = 2 one arrives to a Taub-NUT-type magnetostatic 
solution. This solution can also be generated, in these coordinates, using a well-known 
theorem proposed by Bonnor (see Ref. [42]) from the Taub-NUT vacuum solution. Note 
that for 6 = 0, it reduces to the Darmois metric [1]. 

From the above expressions we can compute the physical quantities associated with 
the discs. We obtain 

e = g - - p^^ (86) 

=Po[l + fcV]^(^'-i), (88) 
ikqipo , . 

" = - 

4a 

^° " {a + l){a + k)^ -{a~l){a-k)^ ^ ' 



where 



is a constant and 

^ = ^ [(« + + ky + {a-l){a- ky] Pr. (91) 

Again we must choose 7 > so that g > and the solution satisfies the strong energy 
condition, po is also a positive constant for these values of 7. 

In Figs. M^a) - ^c) the plots of the quantities e, p^,, Pr, and a are presented for 
discs with a = 3, 7 = 1.2, and qi = q2 = 0, 0.5, 0.8, and 1.5, as functions of r. As in the 
previous case, we see that the energy density initially is a positive quantity and then 
becomes negative in violation of the weak energy condition. Therefore only the central 
region of these discs has a physically reazonable behavior. The another functions also 
have a similar behavior to the precedent case. We also study Zipoy-Voorhees-type discs 
for other values of the parameters, but in all the cases we found a similar behavior. 
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In the same way, the relevant quantities of the CRM are shown in the following 
figures for the same values of the parameters, also as functions of r. We also consider the 
case when the two fluids move with equal and opposite tangential velocities. Here the 
tangential velocity v"^ (Fig. W^d)) is always less than the light velocity. One also finds 
that the inclusion of electromagnetic field makes less relativistic these discs. We also 
obtain h'^ (Fig. EJ^e)) always as an increasing monotonous function of r that corresponds 
to a stable CRM for the discs. Finally, in Fig. ^f) we have plotted the energy and 
electric charge densities e± and a± in the particular case when the two streams have 
equal pressure or equal energy density. We see that a± has a similar behavior to a, 
and e± to the energy density e. Therefore, for Zipoy-Voorhees-type fields we can build 
counterrotating thin disc sources only with a well-behaved central region. 



4-.3. Discs from a Bonnor-Sackfield-type solution 



The third family of solutions considered is a Bonnor-Sackfield-type solution which is 
given by 

2 (92) 



A 

At 



(a + l)e'^=°* - (a - l)e-T™* 
-(7V2)ln \{u^ + l)/{u^ + w'')], 



V2qi\ 



^7 cot u 



~7 cot 



(a + IjeT^^t " 
A^ = \/2k'yq2W, 



(a — l)e 



-7 cot It ' 



(93) 
(94) 
(95) 



where = 1 + 6^, with = Q'f + and 7 is a real constant. Here gi and are again 
the electric and magnetic parameters, respectively, u and w are the oblate spheroidal 
coordinates, related to the Weyl coordinates by 



kuw, 



(96) 



where < m < 00 and < w < 1. In this case we must choose < a < 1 and the 



shell would be located at u = a/k > and w = ^/l~—r^, with k = Vl — As in the 
precedent cases this solution can be generated, in these coordinates, following Ernst's 
method ^\ from the Bonnor-Sackfield vacuum solution [2]. 

From the above expressions we can compute the physical quantities associated with 
the discs. We obtain 



e 

Pr 



Q-Pr- Pip, 
p„[l_A;V]-(7^+l)/2 



1 + 

1 — fc^r^ 
ikqiPo 



Pr 



Pr 



where 



4a 



(97) 
(98) 

(99) 
(100) 
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is a constant and 

ikpo 



:a + l)e^"°*' " + (a - 1)6"''^°*" 



Pr. (102) 



2a2 

Here we also must choose 7 > so that ^ > and the solution satisfies the strong 
energy condition, po is also a positive constant for these values of 7. 

In Figs. [31(a) - ^c) the plots of the above quantities are presented for discs with 
a = 0.5, 7 = 1.3, and qi = q2 = 0, 0.5, 0.8, and 1.5, as functions of r. As in the previous 
cases, we see that the energy density initially is a positive quantity and then becomes 
negative in violation of the weak energy condition. Therefore only the central region 
of these discs has a physically reazonable behavior. The another functions also have 
a similar behavior to the earlier cases. We also study Bonnor-Sackfield-type discs for 
other values of the parameters, but in all the cases we found a similar behavior. 

Again, the relevant quantities of the CRM are shown in the following figures, also 
as functions of r. We also consider the case when the two fiuids move with equal and 
opposite tangential velocities. Here the tangential velocity (Fig. EJ^rf)) is not always 
less than the light velocity. However the inclusion of electromagnetic field can make 
the velocities of the particles smaller than the light velocity. In some cases h'^ (Fig- 
Et^e)) presents strong changes in the slope at certain value of r, which means that there 
is a strong instability there. For other values of the parameters h"^ is an increasing 
monotonous function of r that corresponds to a stable CRM for the discs. 

Finally, in Fig. [3](/) we have plotted the energy and electric charge densities e± 
and a± in the particular case when the two streams have equal pressure or equal energy 
density. We see that e± has a similar behavior to the energy density e. Therefore, for 
Bonnor-Sackfield-type fields we can build counterrotating thin disc sources only with a 
well-behaved central region. 

4.4- Discs from a charged and magnetized Darmois solution 

Finally, other family of solutions to the Einstein-Maxwell equations is a charged and 
magnetized Darmois solution, which is given by 

a^x^ — — 1 , , 

[ax + 1)^ — o^y^ 



A = 21n 



ax — by — 1 
a^(x^ — y 



2\ 



(104) 



2V2qiy 

^* = 7 — — (-10^) 

[ax + 1)"^ — b'^y'^ 

V2kq,{l-y'){ax + l) 
a(a2x2 - 6V - 1) ' ^ ' 

where = 1 + 6^, with b"^ = qf + q2- Here qi and q2 are also the electric and magnetic 
parameters, respectively, x and y are the prolate spheroidal coordinates, given by 
(!85|) . Again we must choose a > 1 and the shell would be located at x = a/k > 1 
and y = \/l — r^, with k = ^ oP- — \. For 6 = this solution also reduces to the 
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Darmois metric. Moreover, for g2 = we have a Kerr-type electrostatic solution. This 
solution can be generated, in these coordinates, using a well-known theorem proposed 
by Bonnor (see Ref. |12]) from the Kerr vacuum solution. And for gi = one obtains 
its magnetostatic equivalent. This solution is asymptotically flat and was firts studied 
by Bonnor [13] in a different system of coordinates and describes the field of a massive 
magnetic dipole. 

The physical quantities associated with the discs can now be written as 



Pr 



a 



Pr 



g-Pr- Py,, 

;i + fc2r2)(a2 + FfeV) 

2aa\l + 

(a2 + k%'^r'^)[{aa + kf - kW{l - r"^)] ' 
4y2agi(aa + A:)(l -r2)i/2 



ai^a^ + k'^h'^r'^)[[aa + kY 
\/2q2kr'^ 



k%^{l-r^)] 



Pr 



4a 



where 



4ak 



Q 



a{a'^ + k'^lP'r'^) 



[aa + kY + k%'^{l - 



Pr 



(107) 
(108) 

(109) 

(110) 

(111) 

(112) 



[aa + kf -k%'^{l-r'^) 

In Figs. 111(a) -Hl^c) the plots of the above quantities are shown for discs with a = 1.5 
and qi = q2 = 0, 0.5, 1.0, and 1.5, as functions of r. We see that these discs also have a 
central region well behaved which satisfies the weak and strong energy conditions and 
border region where e < 0, in violation of the weak energy condition. We also computed 
these functions for other values of the parameters, but in all the cases we found a similar 
behavior. 

Equally, the relevant quantities of the CRM are shown in the following figures, also 
as functions of r. We also consider the case when the two fluids move with equal and 
opposite tangential velocities. We see that the inclusion of electromagnetic field can 
make the velocities of the particles smaller than the light velocity. Here [Fig. W^d)] 
can also take negative values. In some cases h"^ [Fig. Hl^e))] presents strong changes in 
the slope at certain values of r, which means that we have a strong instability there, 
and also we find regions with negative slope, showing that the CRM can not be applied 
for these values of the parameters. However for qi = q2 = 0.7 h"^ is always an increasing 
monotonous function of r that corresponds to a stable CRM for this disc. 

Finally, in Fig. S)^/) we have plotted e±, (T+, and (T_ in the particular case when the 
two streams have equal pressure or equal energy density, for values of the parameters for 
which f 2 < 1. We see that e± has a similar behavior to the energy density e. Therefore, 
for this solution we can build counterrotating thin disc sources only with a well-behaved 
central region. 
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5. Discussion 

A detailed study was presented of the counterrotating model for generic electrovacuum 
static axially symmetric relativistic thin discs with nonzero radial pressure. A general 
constraint over the counterrotating tangential velocities was found, needed to cast the 
surface energy-momentum tensor of the disc in such a way that it can be interpreted as 
the superposition of two counterrotating charged perfect fluids. The constraint found 
is completely equivalent to the necessary and sufficient condition obtained in Ref. |31j . 
We show that, in general, it is not possible to take the two counterrotating fluids 
as circulating along electrogeodesics neither take the two counterrotating tangential 
velocities as equal and opposite. We also have obtained explicit expressions for the 
energy densities, current densities, and velocities of the counterrotating streams in terms 
of the energy density, azimuthal pressure, and planar current density of the discs, that 
are also equivalent to the correspondig expressions in Ref. [3T] . 

Four families of models of counterrotating charged discs were considered in the 
present work based on simple solutions to the vacumm Einsteins-Maxwell equations in 
the static axisymmetric case generated by conventional solution-generating tecniques 
[1]. We saw that the presence of an electromagnetic field can make the counterrotating 
tangential velocities of the particles smaller than the light velocity, and also can stabilizes 
againts radial perturbations the CRM. The discs constructed have a well-behaved central 
region which satisfies the weak and strong energy conditions and border region where 
the energy density is negative, in violation of the weak energy condition. 

Finally, the generalization of the counterrotating model presented here to the case 
of rotating thin discs with or without radial pressure in presence of electromagnetic 
fields is being considered. 
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e p^, Pr 




(e) (/) 

Figure 1. For Chazy-Curzon-type fields we plot, as functions of r, (a) the energy 
density e and (b) the pressures pip (upper curves) and pr for discs with 7 = 2 and 
9i = 92 = (top curves), 0.5, 1.0, and 1.5 (bottom curves), (c) the electric charge 
density a for discs with 7 = 2 and qi = q2 = (axis r), 0.2, 0.3, 0.5, and 1.5 (top 
curve), (rf) v"^ for discs with 7 = 0.7 and qi = q2 = 2 (top curve), 3, 5, and 10 (bottom 
curve), (e) the angular momentum for discs with 7 = 0.7 and gi = 52 = 2 , 3 (sharp 
curves), 5, and 10 (top curve), and (/) e± (lower curves) and a± for discs with 7 = 0.7 
and qi = q2 = 5 and 10 (top and bottom curves, respectively). 
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ic) 



(d) 
e±, -a± 



(e) (/) 

Figure 2. For Zipoy-Voorhees-type fields we plot, as functions of r, (a) the energy 

density e and (6) the pressures p^p (upper curves moved upwards a factor of 0.5) and pr 
for discs with a = 3, 7 = 1.2 and qi = q2 = (top curves), 0.5, 0.8, and 1.5 (bottom 
curves), (c) the electric charge density a for discs with a = 3, 7 = 1.2 and gi = (72 = 
(axis r), 0.5, 0.8, and 1.5 (top curve), (d) for discs with qi = qi = ^ (top curve), 
1, 2, and 5 (bottom curve), (e) the angular momentum /i^, and (/) e± for discs with 
a = 3, 7 = 1.2, and gi = (j'2 = (bottom curves), 0.5, 0.8, and 1.5 (top curves), and 
o'i (upper curves scaled by a factor of 10) for qx = <li=^ (axis r), 0.5, 0.8, and 1.5 
(top curve) and the same values de a. and 7. 
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Figure 3. For Bonnor-Sackfield-type fields we plot, as functions of r, (a) the energy 

density e and (6) tlic pressures ptp (upper curves) and pr for discs with a = 0.5, 7 = 1.3 
and 5i = 52 = (top curves), 0.5, 0.8, and 1.5 (bottom curves), (c) the electric charge 
density a for discs with a = 0.5, 7 = 1.3 and qi = q-z = (axis r), 0.5, 0.8, and 1.5 
(top curve), (d) for discs with a = 0.9, 7=1 and (ji = 92 = 3 (bottom ciuvc), 5, 8, 
and 10 (top curve), (e) the angular momentum for discs with a = 0.9, 7 = 1, and 
?i = 52 = 3, 5 (sharp curves), 8, and 10 (bottom curve), and (/) e± (lower curves) 
and <j± for discs with a = 0.9, 7=1, and qi = q2 = 8 and 10 (top and bottom curves, 
respectively). 
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e p^, Pr 




Figure 4. For charged and magnetized Darmois fields we plot, as functions of r, (a) 

the energy density e for discs with a = 1.5 and (71 = 92 = (curve with crosses), 0.5, 
1.0, and 1.5 (bottom curve), (6) the pressures (upper curves) and Pr for discs with 
a = 1.5 and qi = q2 = 0, 0.5, 1.0, and 1.5 (top curves at the border of the discs), (c) 
a (upper curves) and j for discs with a = 1.5 and gi = (72 = (axis r), 0.5, 1.0, and 
1.5 (top curves), (b) for discs with a = 1.5 and qi = q2 = (top curve), 0.5, 0.6, 
0.7, 0.8, and 1.0 (bottom curve), (e) the angular momentum for discs with a = 1.5 
and qi = q2 = 0, 0.5, 0.6 (sharp curves), 0.7 (top curve), 0.8, and 1.0 (bottom curve), 
where the last curves have been scaled by a factor of 20, and (/) e±, c7+ (upper curves), 
and a_ (lower curves) for discs with a = 1.5 and qi = q2 = 0.7 and 0.8 (top, bottom 
and bottom curves, respectively). 
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